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ABSTBACT 

This  thasis  utilizes  historical  wearmetal  data  from  the  Ft. 
Ord  Spectrcmetric  Oil  Analysis  laboratory  to  propose 
wearmetal  level  and  trend  guidelines  for  equipment  powered 
by  the  Continental  LD/LDS/LDT  465/465-1  engine.  The  msthod- 
ology  proposed  for  determining  trend  guidelines  requires 
data  giving  the  parts-per- million  (PPM)  level  of  a  given 
wearmetal   at  a    known    time   since  oil   change. 
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I-    IHTROPgCTIOg 

During  the  operaticr.  of  any  place  of  mechanical  equip- 
ment, a  certain  amount  of  abrasion  occurs  between  metal 
parts  cciing  in  contact  with  one  another.  The  result  of 
this  abrasion  is  the  deposition  of  minute  quantities  of 
metal    (wearmetals)     in   the  lubricating    oil   of   the   equipment. 

It  was  recognized  by  the  railroads  as  far  back  as  the 
1940* s  that  monitoring  the  concentration  cf  wearmetais  in 
the  lubricating  oil  could  be  useful  in  determining  the 
mechanical  condition  of  diesel  engines  and  forewarn  the 
breakdown   of  engine    components. 

Later  in  the  mid  1950* s,  the  Navy  began  a  trial  moni- 
toring program  at  the  Naval  Air  Rework  Facility  in  Pensacola 
to  determine  if  oil  analysis  techniques  could  be  useful  for 
detecting  th.9  abnormal  operation  of  aircraft  engines. 
Following  the  success  of  the  Navy  effort,  the  Army  first 
applied  the  techniques  to  the  monitoring  of  reciprocating 
engine  helicopters  in  1959.  An  evaluation  of  oil  analysis 
programs  was  initiated  by  the  Air  Force  in  1962  and  resulted 
in  the  establishment  of  an  Air  Force  oil  analysis  testing 
program    in    1964. 

Consolidation  of  Department  of  Defense  oil  analysis 
programs  was  initiated  in  1969  with  the  DOD  Equipment  Oil 
Analysis  Program  (SOAP)  and  later  the  Joint  Oil  Analysis 
Program  (JOAP)  .  The  purposes  of  this  consolidation  were  to 
effect  uniformity  and  standardization  in  oil  analysis 
program  sguipment,  standards  and  techniques,  consolidation 
cf  procurement  requirements  for  oil  analysis  equipment  and 
the  centralization  cf  responsibility  for  technical  manage- 
ment   and  overall   program   surveillance. 


The  basic  assumptions  of  ths  spactrometric  analysis  of 
wearmetals  is  that  during  the  normal  operation  of  a  piece  of 
equipment ,  the  interacting  of  parts  results  in  the  produc- 
tion of  minute  quantities  of  metals.  This  wearraetal 
production  cccurs  at  a  constant  rate  (as  a  function  cf  time) 
for  properly  operating  pieces  of  equipment.  However,  as  a 
part  enters  into  its  failure  mode,  a  much  higher  level  of 
wearmetal  particles  may  be  deposited  into  the  lubricating 
oil.  An  illustration  of  this  process  for  a  hypothetical 
piece    of    equipment   is   shown    in   Figure    1.1. 


EEM 


Break 


Abnormal 


HOUES 


Figure   1.1        Hypothetical   Wearmetal   Production, 


There  are  limitations  to  the  types  of  failure  which  can 
be  identified  through  spect rcmetric  analysis.  The  types  cf 
failure  which  can  be  detected  are  failures  which  are  accom- 
panied by  wearmetal  production  and  which  occur  at  a  slow 
enough  rate  to  allow  maintenance  action  after  identification 
of  a  possible  failure.  Some  of  the  types  of  unidentifiable 
failure   are: 


1)       failures    which   produce  large  wearmetal   particles    visi- 
ble  tc  the   naked  eye 


2)  failures    which    cccur  too   rapidly  to    be   detected   by   the 
current  oil  analysis  procedures 

3)  failures    which    cccur  without  appreciable   production   of 
wearmetals 

To  identify  impending  failure  conditions,  wearmetal 
levels  in  an  oil  sample  from  the  equipment  studied  are 
determined  through  spectrom etric  analysis.  Next,  the  labora- 
tory personnel  calculate  a  10  hour  PPM  increase  figura  by 
finding  the  change  in  PPM  between  samples,  dividing  that 
figure  by  the  hours  between  samples  (giving  a  rate  of 
increase  per  hour)  and  multiplying  by  10  (giving  the  rate  of 
increase  for  10  hours).  A  guideline  chart  giving  wearmetal 
PPM  criteria  for  normal,  marginal,  high  and  abnormal  ranges 
and  fcr  abnormal  10  hour  trends  is  consulted  for  the  appli- 
cable type  of  equipment.  These  guidelines  are  not 
recommended  wearmetal  values,  but  indicate  ranges  of  values 
within  which  the  EPM  levels  of  similar  equipment  have 
historically  fallen. 

A  determination  is  made  as  to  whether  the  amount  of 
wearmetals  in  the  sample  is  considered  in  the  normal, 
marginal,  high,  or  abnormal  range  and  whether  the  PPM 
increase    in    the      10    hcur   period   is   abnormal.  Based  on   the 

range  of  the  sample,  the  range  of  the  previous  sample,  the 
trend  reading  and  whether  the  oil  sample  was  a  routine 
submission  or  the  result  of  a  laboratory  request,  a  recom- 
mendation is  made  utilizing  the  Decision  Making  Guidance 
Table  (Table  6-1)  of  TM  38-301.  Basically  the  table  calls 
for  increased  sampling  if  there  is  an  indication  of  abnormal 
operations.  Continued  abnormal  readings  result  in  the  dead- 
lining  of  the  equipment  and  a  recommendation  that 
maintenance  action  be  initiated  to  determine  the  cause  of 
the    high    level    of   wearmetals.  In   some    cases  a    combination 

of    abncrmal    wearmetal     readings  can   be    used      to    pinpcint   the 


specific      area    of      the  equipment      which   is      probably    in      the 
process   cf    failing. 

The  maintenance  personnel  must  then  provide  followup  to 
the  oil  analysis  laboratory  stating  whether  -heir  indication 
of  impending  failure  was  correct.  (known  as  a  'hit1)  or 
incorrect    (known  as    a    'miss'). 

Thus,  through  the  combined  efforts  of  the  oil  analysis 
laboratory  and  the  maintenance  personnel,  many  impending 
failures  can  be  detected  and  corrected,  preventing  further 
damage  (cr  perhaps  destruction)  of  the  equipment  and  serious 
injury   to   -he  personnel   operating  it. 
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II.    NAT ORE    OF    THE   PROBLEM 

Extension  of  the  oil  analysis  program  to  Army  non- 
aeronautical  equipment  is  fairly  recent,  therefore  there  are 
a  number  of  engines  in  the  Army  inventory  which  are  still 
lacking  wearmetal  and  trend  guidelines.  The  DARCOM  Materiel 
Readiness  Support  Activity  at  Lexington,  Ky.  provided  a  list 
of  engines,  used  in  Army  ground  equipment,  which  were 
lacking  published  guidelines.  Of  the  engines  listed,  the 
one  which  was  most  widely  used  (in  terms  of  types  of  equip- 
ment of  different  nomenclature)  was  the  continental 
LD/LBS/LCT  465/465-1  (henceforth  referred  to  as  the 
LDS-465).  This  is  a  six  cylinder  diesel  engine  which  pcwers 
a  large  number  of  the  Army»s  2  1/2  Ton  and  5  Ton  trucks.  A 
list  of  the  types  of  vehicles  powered  by  the  engine  is  at 
Appendix   A. 

The  goal  of  the  analysis  was  to  examine  historical  data 
(from  DD  Form  2027-  Cil  Analysis  Record)  on  a  large  sample 
of  different  vehicles  operated  in  different  parts  of  the 
country,  which  were  equipped  with  the  LDS-465  engine.  From 
this  analysis  it  was  hoped  to  derive  guidelines  for  what 
could  be  ccnsidered  normal,  marginal,  high  and  abnormal 
concentrations  of  wearmetals  in  the  lubricating  oil.  These 
guidelines  were  to  be  general  (i.e.  applicable  to  all  types 
of  equipment  using  the  subject  engine)  in  keeping  with  the 
format  of  the  guidelines  published  in  TM  38-301  and  TM 
38-30  1-1.  In  addition,  it  was  desired  to  determine  the 
relationship  between  operating  hours  of  the  piece  of  equip- 
ment and  the  rate  of  wearmetal  production.  This  information 
was  to  be  used  to  develop  guidelines  on  what  was  to  be 
considered  an  abnormal  increase  in  wearmetal  PPM  over  a  10 
hour    operating    period. 
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Due  tc  time  ccnstiaints  in  conducting  the  analysis,  Army 
wide  data  could  not  be  obtained.  The  scope  of  the  analysis 
was  therefore  limited  to  a  sample  of  data  on  equipment 
serviced  by  the  Ft.  Ord  Spectrometric  Oil  Analysis 
Laboratory.  This  limitation  on  the  scope  of  the  analysis 
also  placed  a  limitation  on  the  applicability  of  the  results 
in  that  they  are  only  representative  of  the  equipment  at  Ft. 
Ord  and  are  based  on  a  relatively  small  sample,  thus  exten- 
sion cf  the  results  to  Army  wide  equipment  shculd  be 
avoided.  The    methodology,        however,         may      be    useful     in 

analyzing   a     much   larger      set    of      data   in      order   to      provide 
guidelines   applicable   to   Army    wide   equipment. 
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III.    ANALYSIS    METHODOLOGY 

A.       DATA    COLLECTION 

Historical  oil  analysis  records  {DD  Form  2027)  for  the 
period  January  1980  through  January  1983  were  examined  for 
76  vehicles  from  Active  Army  units  located  at  Ft.  Ord.  The 
records  studied  were  for  2  1/2  Ton  Cargo  trucks  of  the  type 
M35A2  or  M35A2C.  A  list  of  serial  numbers  for  the  vehicles 
examined   is    at    Appendix      3. 

The  M35A2  and  M35A2C  are  tactical  vehicles  and  as  such, 
the  guidance  on  routine  sampling  intervals  is  once  every  two 
months.  Samples  are  taken  whils  the  oil  is  warm  either  by 
the      use    cf     a    dip      tube   or      by   a      drain   outlet.  Detailed 

descriptions  of    the    sampling   mathods      are    found    in   TM    38-301 
or    TB    43-0211. 

An  Oil  Analysis  Bequest  (DD  Form  2026)  is  filled  cut 
giving  such  information  as  equipment  identification  data, 
time  oil  sarcple  taken  and  method  used  (drain  or  dip  tube), 
hours  since  oil  change,  hours  since  overhaul,  and  any  state- 
ments about  unusual  operation  of  the  equipment.  The  samples 
are  then  sent  to  the  Oil  Analysis  Laboratory  at  Ft.  Ord 
where  a  spectrcmetric  analysis  is  performed  to  determine 
wearmetal    levels   in    the    sample. 

The  analysis  is  accomplished  by  the  use  of  the  Baird 
Emission  Spectrometer  Model  A/E350-3.  A  small  portion  of 
the  oil  sample  is  placed  in  the  spectrometer.  A  carbon  disk 
electrode  rotates  through  the  sample,  picking  up  a  film  of 
oil  which  is  burned  by  an  electric  arc  between  a  fixed  elec- 
trode and  the  rotating  disk.  The  energy  radiated  by  the 
burning   oil   is    separated   into      its   component   wavelengths   and 


13 


the  intensity  measured.  This  is  compared  to  the  intensity 
measured  frcm  calibration  samples  to  determine  the  parts  per 
million  (PPM)  levels  of  the  various  wearmetais  in  the 
sample.  There  is  a  degree  of  variablility  associated  with 
the  spectrometer  reading.  TM  38-301  states  the  spectrometer 
has  a  tolerance  of  ±1.0  PPM;  however  it  seems  more  likely 
the  variability  is  a  function  of  the  magnitude  of  the 
reading   as   suggested    by  the   spectrometer   manufacturer. 

After    spectrometric  analysis,        the   following    information 
is   then   entered    into   the   DD      Form   2027    (Oil    Analysis    Record) 
for   the   sampled    piece  of   equipment. 

1.  Sample  number     (laboratory   assigned) 

2.  Data   index    (    data   correction  or    maintenance    feedback) 

3.  Julian  date   saiple    analyzed 

4.  Laboratory  response   time 

5.  Last   laboratory  recommendation 

6.  Hcuis    since  overhaul 

7.  Hours    since   oil  change 

8.  Season   for   the   sample 

9.  Wearmetal    levels 

10.  Post   analysis    data    (    lab   recommendations  and    mainte- 
nance   follcwup) 

None  of  the  records  studied  indicated  laboratory  recom- 
mendations for  maintenance  action  and  thus  no  reports  of 
maintenance     followup     existed.  It    was      also      found     that 

approximately  25^  of  the  historical  records  lacked  data  for 
hours      since  oil      change   and      hours      since   overhaul.  This 

precluded   any      systematic  analysis  of      the   entire      sample   of 
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wearmetal   levels     as   a  function      of   operating   hours.  As   a 

result,      to   determine  wearmetal   guidelines,      an   analysis    was 
conducted    which    did    net   require   operating    hour   information. 

The  suggested  oil  sampling  intervals  were  also  not 
adhered  tc  in  all  cases.  As  a  result,  the  number  of  records 
available  fcr  each  equipment  during  the  time  period  studied, 
ranged   frcm    14    to    17    samples. 

B.       DETEBBINATION   OF    HBARMBTAL    LEVEL    GUIDELINES 

Faced  with  these  constraints,  the  method  of  analysis 
chosen  was  to  look  at  all  of  the  PPM  levels  available  for 
the  76  vehicles,  plot  a  frequency  distribution  and  examine 
the    quantiles  from    the  empirical   distribution. 

The  director  of  the  Ft.  Ord  Oil  Laboratory  provided 
guidance  on  the  major  wearmetals  of  concern  for  the  LDS-465 
engine.      They   were: 

1.  Iron    (Fe) 

2.  AluiEir.um     (Al) 

3.  Chrcmium    (Cr) 

4.  Copper    (Cu) 

5.  Lead    (Pb) 

Thus  the  scope  of  the  study  was  limited  to  these  five 
wearmetals.  The  wearraetal  data  was  read  by  the  FORTBAN 
program  WEARMTL  (see  Appendix  F  for  computer  program  list- 
ings) which  produced  an  output  frequency  histogram  for  each 
wearmetal  alcng  with  a  number  of  sample  statistics.  In  all, 
a  total      cf    1255   sample      points  were    available.  The    60th, 

70th,    80th,    90th,    95th,    96th,    97th,    98th,    and   99th    quantiles 
were    determined    from      the  data.        No    decisions      were   made  on 
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what  labels  (i.e.  normal,  marginal,  high,  or  abnormal)  to 
place  on  these  guantiles,  but  rather  to  provide  this  infor- 
mation tc  the  Oil  Analysis  Laboratory  director  for  use  in 
his      own   subjective      judgement.  The      guantiles    found      are 

listed    in   Table    I.    Output  histograms    are   at   Appendix   C. 


■ 

TABLE    I 

Sel<= 

jcted 

learmetal 
Quan tile 

Quantiles 

Wearmetal 

Iron 

Alumiui 

Chromium 

Ccpper 

Lead 

6  0th 
94 
20 
9 
16 
39 

70th 
108 
25 
10 
19 
48 

80th    90th 
124      149 
30         38 
12         15 
24         31 
61         76 

95th    96th 
175       180 
46         47 
18         18 
38         41 
91         95 

97  th 

200 

50 

19 

46 

104 

98th 

215 

53 

21 

51 

1  16 

99  th 

238 

57 

24 

73 

137 
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C.       DETERMINATION    OF    10    HOO B    TREND    GUIDELINES 

In  addition  to  wearmetal  levels,  a  major  concern  in  the 
oil  analysis  program  is  wearmetal  trends.  Trend  guidelines 
are  established  for  a  standard  10  hour  period  and  give  what 
is  considered  to  be  an  abnormal  increase  in  wearmetal  levels 
for      that   period.  In     the   study      performed      by    the      ARINC 

Corporation  for  the  0.  S.  Air  Force  [ Ref .  1,  pp  10-14]  the 
method  of  piecewise  or  segmented  regression  was  used  to 
determine   the   point      of   onset   of    failure.  In    this    method, 

wearmetal  concentration  data  on  a  number  of  similar  pieces 
of  eguipient  is  divided  into  two  groups,  that  data  repre- 
senting     normal    wear      and  that      representing   abnormal      wear. 
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The  hours  associated  with  each  sample  are  converted  from 
hours  since  oil  change  to  hours  prior  to  detection  of 
failure.  This  indicates  the  necessity  of  maintenance  feed- 
back to  determine  or  verify  a  point  of  failure  detection. 
Straight  lines  are  then  fit  through  both  sets  of  data  using 
the    procedure  detailed     by    Hudson   [Ref.       2].  The    point   of 

intersection  of  the  two  lines  is  taken  to  be  the  time  prior 
to  detection  at  which  failure  commenced  (see  Figure  3.1). 
The  slope  of  the  line  following  this  point  could  then  be 
used    to   calculate  abnormal    10    hour  trends. 


PEH 


Normal 


B (Normal) 


Abnormal 


B  (failure) 


Failure   commences 


Operating  hours   prior  to    dstection 


Figure   3.  1        Fitting  Data   By  Piecawise  Regression, 


Since  no  maintenance  feedback  was  available,  a  technique 
was  proposed  which  would  indicate  abnormally  high  metal 
trends  without  requiring  maintenance  information  on  actual 
failure   conditions.      This   method   assumes   that: 

1.  Wearmetal    levels  may    depend   on    operating   hours   of  the 
equipment. 

2.  A    possibly   different    relationship   exists  for    the   nor- 
mal   and  abnormal  periods    of  operation. 
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3.      The   rate    of  wearmetal    production   is    higher   during 
failure  than    during  normal    operation. 

Instead  of  fitting  piecewise  regression  lines  to  the 
data,  the  proposed  methodology  treats  the  slope  of  each  line 
segment  (whose  end  points  are  formed  by  successive  wearmetal 
readings)  as  an  estimate  of  the  rata  of  wearmetal  produc- 
tion. An  illustration  of  this  technique  is  given  in  Figure 
3.2.         Each    dotted   line   segment      provides   an    estimate   of    the 


HCDRS    SINCE    OIL    CHANGE 


Figure    3.2        Piecewise  Estimation   10    Hour    aearaetal   Trends. 


production      rate.  If      the    wearmetal      production      rate     is 

linear,    the   relationship   is    of   the  form: 


Where  : 


Y  (i)        =    B(0)    +    B(1)*X(i)    ♦    S(i) 

X  (i)  =  hours  since   oil    change    for   the    ith   sample 

Y  (i)  =   wearmetal    reading  for   the    ith   sample 
B  (0)  =   intercept    of   the  regression   line 
B(1)  =   slope   of  the   regression   line 
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E  (i)  =  deviation  of  the  actual  wearmetal  reading 
from  the  regression  line. 

The   above  definitions  apply   throughout   the  rest   of   the 
thesis. 

The  slope  of  each  successive  line  segment  is  found  by  the 
formula 

slope  =  (Y(i  +  1)  -Y(i)  )/(X(i+1)-X(i)) 

substituting  the  regression    equations: 

slope    =    £  <0)  +  (3  (1)  *X  (i+1))  +E(i»1)  -B(0)-  (B  (1)*X(i)  )-E  (i) 

X(i  +  1)  -X(i) 

=  £(1)    +    (E(i  +  1)-E(i))  /(X(i  +  1)-X  (i)) 

Sines  the  regression  slope  B(1)  is  constant  for  any  given 
set  of  data,  each  pairwise  slope  will  deviate  from  the 
actual   slope  by    a   factor 

(E  (i+i)-E(i))/(X(i+l)  -x<i)) 

Thus,  the  •goodness*  of  this  method  is  determined  by  the 
magnitude  of  the  variance  in  the  residuals  and  the  length  of 
time  between  wearmetal  samples.  Samples  with  small  vari- 
ances in  their  error  terms  and  large  sampling  intervals  will 
approximate  the  regression  slope  much  better  than  those  with 
large   variances    and    small  sampling  intervals. 

The  next  step  in  the  methodology  is  to  convert  each  of 
these  slopes  to  a  PPM  increase  for  a  10  hour  period.  This 
is   done    ty   multiplying  the    slope    by    10: 

Increase    in   PPM    =    1  0*  (B(  1)  +  ( ( (E  (i  +  1)  -  E  (i)  )  /  (X  (i+1)  -X  (i)  ) )  ) ) 

in    10    hours 

Thus  each  of  the  slope  segments  (regardless  of  the  time 
interval    between  the    samples  )       is      converted   to    a   change   in 
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PPM  which  would  be  expected  during  any  10  hour  period  within 
this  interval.  The  10  hour  PPM  rate  of  change  measurements 
are  then  grouped  together  in  a  frequency  histogram.  If  the 
rate  of  wearmetal  production  is  higher  during  failure  than 
during  ncrmal  operation,  then  the  upper  quantiles  of  the  FPM 
rate  cf  change  distribution  should  represent  these  samples 
from    the    failure   mode  of   the   equipment. 

This  methodology  was  applied  to  the  Fort  Ord  data  (for 
those  data  elements  which  had  hours  since  oil  change  infor- 
mation available)  through  the  use  of  the  FORTRAN  program 
ANALYSIS  (Appendix  F) .  The  logic  used  in  computing  the 
slopes    was: 

1.  If   no   hours  since   oil    change  data    was   available,    the 
sample  was   discarded. 

2.  If   hours    since    oil   change    was   greater   than    the   pre- 
vious  reading,    the   slope   was  determined   by: 

(Y(i  +  1)-Y(i))/(X(l+1)  -X(i)  ) 

3.  If   hours    since   oil   change   was   less   than   that   of   the 
previous    sample,   it   was    assumed    an   oil   change    took 
place    and    thus    the   slope    was  calculated    by: 

Y(i)/X(i) 

Since  approximately  25%  cf  the  records  lacked  time  since 
oil  change  data,  only  8  76  trend  values  were  computed. 
Histograms  cf  t-he  trend  frequencies  are  at  Appendix  D.  The 
results  obtained  shoned  a  large  variance  in  the  PPM  rate  of 
change  levels.  This  may  have  been  due  to  errors  in  data 
reporting  since  a  large  factor  in  the  variance  was  the  exis- 
tance  of  a  few  extremely  large  or  small  measurements.  As  in 
the  case  of  the  wearmetal  levels,  selected  quantiles  were 
obtained  and  are  shown  in  Table  II. 
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TABLE    II 
Selected   10  Hour  Trend   Quan-ciles 

Q  uantile 

Wearmetal 

Ircn 

Aluminum 

Copper 

Chromium 

Lead 

60th    70th    80th    90xh    95th    96th    97th 
6.0       9.3      15.7    31.9   66.9    86.3    140. 
1.2       2.0        3.5      8.5   20.0    23.8    30.0 
1  .1        1.9         3.3      7.5    15.9    20.9    26.9 
0.5       0.9         1.5      3.3      8.0    10.0    13.0 
2.4       4.5         8.8    18.3    37.5    42.4    52.5 

98th 
190. 
41  .3 

35.0 
20.0 
30.0 

99  th 
230. 
70.0 
45.9 
23.3 
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IV.     IlIDSTBATION    OF    TREND    ANALYSIS    TECHNIQUES    OSING 

SIMOIATED    DATA 

The  technique  proposed  for  analyzing  trends  was  not 
adequately  illustrated  by  the  actual  data  for  the  76  vehi- 
cles in  the  sample  due  to  the  suspect  validity  of  the  hours 
since  oil  change.  To  further  illustrate  the  technique,  a 
set  cf  data  was  generated  which  fit  the  assumptions  set 
forth  when  proposing  this  technique.  Those  assumptions 
were: 

1.  PFW   levels   may    te   related   to  the   operating    hours 
since    oil    change. 

2.  The   wearmetal    production   process    may   be   different   for 
the   normal  and    abnormal    modes   of   operation. 

3.  The   rate    of   wearmetal    production   is    higher    during 
failure  than    during  normal    operation 

Thus    simulated   wearmetal    data    was    generated  by    assuming: 

Y  (i)     =    B(0)     ♦    B  (1)  *X(i)     ♦    E  (i) 

Where   the   variables    are   defined   as  in    chapter    III. 

Since  the  proposed  technique  approximating  the  regression 
line  slope  by  the  individual  piecewise  slopes  was  indepen- 
dent of  the  intercept  term  (i.e.  B(0)),  it  was  only 
necessary  tc  derive  estimates  for  B(1)  and  the  variance  of 
the    error   terms. 

To  determine  a  representative  value  to  use  for  B(1), 
records  fcr  the  76  trucks  were  examined  to  determine  whether 
any  of  the  vehicles  had  wearmetal  data  for  iron  which  exhib- 
ited  the    following    characteristics: 
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1.  Continuous,   relatively   uniform    increases   in    hcurs 
since    oil    change  data 

2.  No   apparent   gaps  in  the    data 

3.  No   abnormal  drcps   in    PPM   levels    without   an    accompa- 
nying  record   of  oil  change 

4.  At   least    8   data   points   fitting    the   above  rsquirements 


A  total  of  ten  vehicles  met  the  above  criteria.  The  data 
from  these  vehicles  was  entered  into  x.he  APL  program  REGRESS 
[Ref.  3],  which  fit  a  least  squares  linear  regression  line 
to  the  data.  A  plot  of  the  residuals  was  examined  tc  deter- 
mine if  there  were  any  apparent  violations  of  the  assumption 
of  normality.  In  addition,  the  R2  values  were  examined  to 
determine  if  there  was  a  strong  relationship  between  the 
dependent  (PPM  levels)  and  independent  (hours  since  oil 
change)  variables.  Those  samples  which  had  R2  values  of  .9 
or  greater  and  did  not  viola-e  the  subjective  test  for 
normality   of  residuals  were    chosen.  A   total   of    four    vehi- 

cles   met    these      criteria.         The   iron    3(1)         values    for   xhsse 
vehicles    ar=  shown    in   Table    III. 


TABLE    III 
Regression   Line  Slopes 

Vehicle  number  B{1) 

3903458  .3648 

3818363  .4277 

385C984  .5141 

3903880  .3306 
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This  technique  is  not  meant  to  represent  a  statistically 
valid  method  for  determining  the  behavior  of  the  wearmetal 
data.  It  is  meant  only  to  provide  a  general  idea  of  the 
magnitude  of  the   wearmetal    production    rates. 

The  variance  was  assumed  to  be  1  for  the  sake  of 
simplicity. 

Thus  tc  generate  the  data,  all  that  was  necessary  was  to 
determine  the  time  at  which  the  sample  was  to  be  taken, 
multiply  it  by  the  slope  (B(1))  and  add  on  an  error  term 
which  was  distributed  N (0,  1)  .  For  illustrative  purposes, 
times  were  generated  at  exact  ten  hour  intervals  (i.e.  10, 
20,  30  etc.  hours)  .  The  error  terms  were  generated  using 
the  LIRANEOM  II  [Ref.  4]  random  number  generating  package. 
To  simulate  abnormal  operations,  data  was  generated  which 
had  a  slope  twice  that  of  the  normally  operating  equipment. 
Error  terms  were  generated  in  the  same  manner  as  for  the 
normally  operating  equipment.  The  data  generated  for  the 
four  pieces  of  equipment  were  combined  and  a  frequency 
distribution  (histogram)  was  formed  througn  use  of  the 
FORTRAN  program  ANA1G2N  (Appendix  F)  .  Five  hundred  fifty 
•normal'  data  elements  and  fifty  'abnormal*  data  elements 
were  generated  for  each  of  the  four  values  of  B(1).  The 
output    frequency   histogram    is   at    Appendix    S. 

Upper  quantiles  of  the  distribution  were  determined  as 
was  dene  with  the  wearmetal  levels.  The  quantiles  found  are 
in    Table    IV. 

A  problem  may  arise  any  time  data  from  different  equip- 
ment is  grouped  together.  Unless  the  assumption  that  all 
equipment  of  the  sane  type  operate  in  the  same  manner  is 
true,  one  could  identify  a  measurement  as  abnormal  in  terms 
of   the      entire    population      when   in      fact      it    is      within      the 
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TABLE    IV 
Selected  Quantiles  of    10    Hour  Wearnetal   Production 

QOANTILE 
WEARMETAL  60th   70th   80th    90th    95th    96th    97th    98th      99th 
Ircn  4.69    5.22   5.84    6.97    8.18    8.83    9.46    10.14    10.81 


normal  operating  range  of  the  equipment  from  which  the  meas- 
urement was  taken.  To  preclude  this  happening  it  would  be 
preferable  if  guidelines  could  be  generated  for  each  equip- 
ment based  en  the  operating  history  of  -chat  equipment  alone. 
An  example  of  the  benefits  of  this  method  is  shown  using  the 
generated   data    from    vehicle    3903458. 

Since  the  error  terms  generated  were  distributed  Normal 
(0,1),  the  distribution  of  the  PPM  reading  for  each  gener- 
ated sample  should  have  been  Normal  (10*3(1),  2).  This 
follows  since  the  generated  intervals  were  10  hours  there- 
fore each  10  hour  PPM  reading  would  be  10  times  the  slope  of 
the  regression  line  (which  was  constant)  plus  the  difference 
in  the  pairwise  consecutive  error  terms.  Since  the  mean  of 
the  error  terms  tas  0,  the  mean  of  the  PPM  readings  would  be 
0+10*E(1),  while  the  variance  would  be  the  variance  cf  B(1), 
which  was  zero,  plus  the  variance  of  the  difference  in  the 
error  terms  which  was  2.  Thus  it  would  be  expected  that  the 
frequency  distribution  for  each  generated  sample  would  be  a 
combination  of  two  normal  distributions  offset  from  each 
other  by  10  times  the  difference  in  their  slopes  (means).  A 
histogram  of  the  data  generated  from  the  S(1)  value  for 
vehicle  3903458  (Appendix  E)  shows  this  to  be  the  case  by 
the    distinct  bimodality   of    the   distribution. 
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Thus  if  tha  assumptions  of  the  wearmetal  bahavicr  during 
normal  and  abnormal  operations  hold,  and  the  difference  in 
wearmetal  production  rates  during  normal  and  abnormal  opera- 
tions is  large  in  comparison  -to  tha  variance  in  the  sample 
of  normal  readings,  then  it  should  be  quite  apparent  whan  a 
reading    falls   out  of    the   range   of   normal   variability. 
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7.     CONCLUSIONS    AND    RECOMMENDATIONS 

The  determination  of  quantiles  from  the  empirical 
distributions  of  wearmetal  levels  and  of  wearmetal  trends 
was  done  to  provide  the  director  of  the  Ft.  Ord 
Spectrcmetric  Oil  Analysis  laboratory  with  general  guide- 
lines on  what  ranges  cf  values  he  could  expect  for  wearmetal 
levels  and  trends  for  the  Continental  LD/LDS/LDT  465/465-1 
based  on  historical  data  for  equipment  supported  by  his 
laboratory.  The  presentation  of  results  was  different  than 
that  in  TM  38-301-1  in  that  no  attempt  was  made  to  attach 
labels  (i.a.  normal,  marginal  etc.)  to  these  quantiles. 
This  approach  was  taken  to  provide  mere  flexibility  and 
information  with  which  to  make  subjective  judgements  en  the 
condition  cf  a  piece  of  equipment  since  ultimately  any 
recommendation  made  by  an  oil  analysis  laboratory  is  subjec- 
tive based  not  only  en  published  guidelines  for  a  type  of 
equip  lent,  but  also  a  knowledge  of  the  operational  history 
of   the    individual   equipment    in   question. 

The  fact  that  operational  history  is  an  important  factor 
in  the  recommendations  made  by  the  laboratory  indicates  the 
current  practice  of  grouping  all  equipment  of  the  same  type 
together  and  applying  general  guidelines  to  the  composite 
sample  may  not  in  fact  be  adequate  for  ail  individual  pieces 
of  equipment,  being  too  conservative  for  equipment  which 
normally  produce  wearmetal  at  a  relatively  high  rate  cr  too 
liberal  for  equipment  which  produce  wearmetal  at  a  rela- 
tively lew  rate.  This  problem  is  illustrated  in  Figure  5.1. 
In  this  figure,  samples  from  two  hypothetical  pieces  of 
equipment  are  shown.  A  regression  analysis  is  performed  on 
each    sample     individually   giving     regression   lines      1    and      2 
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Figure    5.1         Regression   of  Combined  Samples. 
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respectively.  These  lines  are  constructed  to  minimize  the 
sum  cf  the  squared  deviations  of  each  sample  point  from  the 
line.  Unless  the  two  regression  lines  are  collinear  (an 
extremely  unlikely  situation)  ,  a  regression  analysis 
performed  on  the  combined  data  would  result  in  a  regression 
line  (line  3)  which,  while  being  the  best  fit  for  the 
combined  data,  would  have  a  larger  residual  sum  cf  squares 
than  either  of  the  individual  sample  regression  lines  indi- 
cating the  fit  of  the  regression  line  for  the  combined  data 
was  not  as  good  a  fit  as  that  of  the  individual  samples. 
Additionally,  the  combined  regression  line  would  overesti- 
mate the  wearmstal  production  rate  of  engine  2  and 
underestimate  the  wearmetal  production  rate  of  engine  1. 
Thus  a  guideline  based  on  the  combined  samples  might  fail  to 
identify  an  impending  failure  in  engine  2  or  errcneously 
signal   an   impending    failure    in   engine    1. 

This  same  problem  applies  to  the  methodology  proposed  of 
estimating  the  regression  line  by  the  individual  10  hour 
trend  readings.  While  each  10  hour  trend  reading  could  be 
considered  a  fairly  good  estimate  of  the  slope  of  the 
regressicn  line  for  that  individual  piece  of  equipment,  it 
would  likely  be  a  pocr  estimate  of  the  slope  of  the  regres- 
sion line  obtained  by  combining  the  data  from  all  pieces  of 
equipment  of  the  same  model.  Again,  the  wearmetal  produc- 
tion characteristics  of  the  individual  equipment  would  be 
masked  by  the  combined  characteristics  of  the  other  pieces 
cf    equipment. 

While  developing  individual  guidelines  on  each  piece  of 
equipment  is  much  more  complicated  than  setting  general 
guidelines  for  all  equipment  of  the  same  type,  these  indi- 
vidual guidelines  wculd  allow  the  laboratory  director  to 
make  a  much  more  accurate  recommendation  after  analyzing  an 
cil   sample. 
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I mplementaticn  of  this  methodology  at  the  laboratory 
level  would  require  computational  capability  (i.e.  a 
computer) .  Sines  the  distribution  of  the  10  hour  PPM  trends 
for  an  individual  piece  of  equipment  is  expected  to  be 
normal,  any  desired  quantiles  of  the  distribution  could  be 
derived  from  the  sample  mean  and  standard  deviation  and  the 
standard  normal  tables.  As  additional  data  elements  were 
addad  (as  a  result  of  additional  samples)  these  statistics 
would  provide  an  increasingly  better  estimation  of  the  true 
wearmetal   trend    characteristics   of  tha   individual    equipment. 

The  gecdness  cf  any  predictive  methodology  is  related  to 
the  quality  of  the  data  on  which  the  prediction  is  made. 
Thus  to  derive  the  maximum  benafits  from  the  oil  analysis 
program,  additional  emphasis  must  be  placed  on  the  proper 
collecting  of  oil  samples  and  reporting  of  equipment  oper- 
ating data.  This  will  require  not  only  increased  command 
emphasis,  but  convincing  maintenance  personnel  of  the 
benefits  of  the  program  in  reduced  maintenance  expense, 
increased   readiness    and   greater   safety. 
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APPENDIX  A 
TYPE  VEHICLES  OSING  THE  CONTINENTAL  LD/LDS/LDT-465  ENGINE 


END  ITEM 

M34A2 

M35A2 

M35A2C 

M36A2 

M45A2 

M46A2 

M49A2C 

H50A2 

M50A3 


NOMENCLATURE 


2  1/2  Ton  Cargo  Truck 


M51A2 

M5  2A2 

M54A2 

M54A2C 

M55A2 

M61A2 

H63A2 

M109A31 

M185A3 

M2U6A2 

2  1/2  Ton  Bolster,  truck: 


2  1/2  Ton  truck  chassis 


2  1/2  Ton  Tank  truck,  fuel 


2  1/2  ion  Tank  truck,  H20 


2  1/2  Ton  Tank  truck,  H20 


5  Ton  Dump  truck 


5  Ton  Tractor  truck 


5  Ton  Cargo  truck 


5  Ton  truck  chassis 


2  1/2  Ton  truck,  tractor,  wrecker 


Inst  Reoair  truck 


5  Ton  truck,  tractor,  wrecker 
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M275A2 
M291A2 
M292A2 
H292A5 
M342A2 


2  1/2  Ton  Tractor,  truck 
5  Ton  truck,  Exp  van 

2  1/2  Ton  truck,  Exp  van 
ii         ii 

2  1/2  Ton  Dump  truck 
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APPENDIX    B 
LIST   OF   SERIAL    NUMBERS    OF   ENGINES    STUDIED 


25880A 

12688 

3903358 

3923514 

3929006 

3914728 

3901793 

3902619 

3765142 

3922795 

3827828 

3891139 

3862237 

391707 

3902420 

3S24675 

3090  8 

3903418 

3923  84  57 

3  925  101 

3903458 

3901270 

8805274 

3650984 

39032  12 

3903350 


381 7608 

3929002 

3923271 

3901901 

3920268 

3921828 

390  1793 

20279 

3806178 

3903772 

10305783 

261  97 

390  3451 

390 1904 

3903137 

3924712 

390  1455 

3922157 

391 2945 

3903884 

3939020 

3903993 

4076463 

3903380 

3816622 


2811 

31244 

3903782 

3902143 

3921134 

3923959 

13778 

6905 

382045 

38176592 

980716 

3818461 

380819 

3902362 

24165 

3902430 

79168 

75145 

3914805 

3047329 

3818363 

3900086 

3925768 

3903421 

3929195 
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APPENDIX    C 
FBEQOENCY    HISTOGRAMS    FOfi    WEAfiflETAL    LEVELS 

This      appendix      contains      histograms      of      the      historical 
wearmetal   level    data    for  the   sample    of   engines   studied. 
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APPENDIX    D 
OUTPUT    FOR    TREHD    AHALTSIS 

The  results  of  this  portion  of  the  analysis  were  highly 
variable.  The  existance  of  a  few  very  large  and  very  small 
values  prevented  the  display  of  the  bulk  of  the  results  in 
detail.  Tc     expand     on   the      display      of      -che      non-extreme 

results,  the  10  largest  and  10  smallest  values  were  trun- 
cated prior  to  producing  the  histograms.  The  values 
truncated   for  each    wearmetal   are    listed   in   Table    V. 


I 

j 

TABLE    V 

_,    ...        .._.,... 

I 

I 
I 

Truncated   aearmetal  Trend 

Values 

I 

IRON 

ALUMINUM 

CHROMIUM 

I 

LOW 

HIGH 

LOW 

HIGH 

LOW 

HIGH 

i 

-680 

22  3.  33 

-280 

64 

-70 

22.22 

I 

-450 

23  0 

-190 

70 

-20 

23.33 

i 

-370 

230 

-60 

70 

-12.31 

27.5 

i 

-135 

.71 

23  2 

-4  0 

70 

-10 

30 

i 

-125. 

.83 

25  0 

-30 

75.71 

-10 

40 

i 

-100 

26  3  .  33 

-3  0 

110 

-10 

40 

i 

-100 

27  0 

-28.57 

110 

-10 

50 

i 

-92.' 

38  3.33 

-25 

130 

-10 

50 

i 

-86.67 

610 

-23.33 

130 

-10 

53.33 

I 

i 

-80 

915 

-22.5 

355 

-6.67 

70 

i 
1 
I 

CCEPER 

LEAD 

i 

LOW 

HIGH 

LOW 

HIGH 

i 

-100 

45.46 

-380 

103.33 

-50 

46 

-260 

110 

i 

-30 

50 

-251. 43 

115 

i 

-21. 

43           62.5 

-120 

117 

i 

-19. 

23           70 

-90 

126.67 

j 

-18. 

33           80 

-83.33 

130 

i 

-12. 

67           83.33 

-8  0 

250 

j 

-12. 

5           110 

-60 

252.86 

1 

-12. 

5           160 

-56.67 

275 

I 

i 

-11. 

52         180 

-4  0 

380 

1 

_    _           j 
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Selected   statistics    on  the    10    hour  trends   are   shewn    in    Table 
VI. 


- 

TABLE    ?I 
Trend    Statistics 

I 

STATISTIC 
MEANS 
STE. DEV 
L.    QRTLS 
MEEIANS 
U.    QRTLS 

Iron 
1  1.88 
62.69 
0.13 
3.49 
1  1.76 

WEARMETAL 
Alum           Chrom 
3.08            2.97 
20.81          12.86 
0.00           0.00 
0.71            0.o1 
2.58            2.44 

Copp 

6".  1  1 
0.00 
0.23 
1.11 

Lead 
5.60 
32.  54 
0.00 
1.  37 
6.  19 

_ 
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APPENDIX    2 
O0TP0T   FOB   TREND    ANALYSIS    USING    SIMULATED    DATA 

This  appendix  contains  histograms  of  the  wearmetai  trend 
analysis  for  the  simulated  wearmexal  data.  The  first  histo- 
gram is  for  the  combined  data  of  four  engines  while  the 
second   is   a   histogram   for   an   individual   engine. 
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APPENDIX    F 
COHPOTEB    PBOGBAHS 


WEARMTL**** 
IN    HEAfiHET 

PLOTS    FREQ 
UER     11    APR 
260)  ,FE{126 
INDX  (1260) 

,J 
NAME    (200) 


C    ****    EEOGFAM 
C    PROGRAM    READS 
C    EQUIEMENT    AND 
C    RICHARE    F.    BA 

REAL    HR(1 
♦  FEM260)  , 

INTEGER    I 

DIMENSION 

J  =  C 

1=0 

DO    15   J=1 

READ 

C  WRIT 

C    READ    IN    THE    W 

10       CONTINUE 

1  =  1+1 

BEAD  (5- 110)     HR(I),FE 
C  WBITEJ6,250)     HR(T),F 


AL    LEVELS    FOR    DIFFERENT 

UENCY    HISTOGRAMS 

1983 

0)  ,AL  (1260)  ,CR(  126  0)  ,CU  (1260)  , 


,75 

(5,100)     NA 
E     (6,270)     N 

EAR  METAL    DA 


ME(J) 
AME  (J) 

TA 


C 

c 

c 
c 
c 
c 
c 


5 
WRIT 


INEX  (I)  =J 
IF     ((HR(I 

CONTINUE 

E    WEARMET 


I^f^i'ttf  ^cfi^^cu^i^PBa] 


))  .GT.-99) 


GO    TO     10 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
GENERATE 


10 

j; 

1  1 

HIST 


CALL    HIST 
WRITE    (6  , 
HI 


CALL 


:st 


WRITE    (6, 

call  hist 
writ: 

CALL 

WRITE    (6, 
CALL    HIST 


HIST 


C 
C 

c 
c 

100 
1  10 
200 
210 
220 
230 
240 
245 
C250 
C270 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

STCP 

FORMAT 

FORMAT 

FORMAT 

FORM  ST 

FORMAT 

FORMAT 

FCRMAT 

FORMAT 

FCRMAT 

FORMAT 

ENE 


6, 
1  1 

10 


AL  LEVELS  TO  APL  DATA  FILE 
,*)   EE 
*)  CU 
*)  CR 

*  AL 
,*)      *E 

OGRAMS    OF    WEARMETAL    LEVELS 

G?     (FE,I,15) 

200) 

GP     (AL,I,15) 

210) 

GP     (CB,I,15) 

220) 

GP     (CU,I,15) 

23  0) 

GP     (EE,I,15) 

240) 

,*)     FE 

*)     AL 

*  CR 
*)     CU 

,*)       EB 


10l 

F10.  1) 


A 

6 
1 
1 

1 
1 
1 
1 
1 

ix;aio 


x,  • 

X,' 

x,« 

x,  • 


HISTOGRA 
HISTOGRA 
HISTOGRA 
HISTOGRA 
HISTOGRA 


X,5F10.5) 
X,6F1C.1) 


OF 
OF 
OF 
OF 
OF 


PPM 
PPM 
PPM 
PPM 
PPM 


LEVELS 
LEVELS 
LEVELS 
LEVELS 
LEVELS 


FOR 

FOR 
FOR 
FOR 
FOR 


ALUMI 


IRON' ) 

INUM') 
CHROMIUM1) 
COPPER') 
LEAD') 
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c 

c 
c 


*****    Pfl 

PROGRAM 

RICHARI 

REAL 

*DELC 

*TRND 

*HR  (1 

*FEJ1 

INTE 

DIME 


DO  5 
READ 
WRIT 
CCNT 


*PB  (I 
CCNT 

CCNT 


C 

C 


10 


C 

C 


15 


49 
50 


*PB  (I 


CGRAM 
ANALY 
I.    BA 

DELH 
R(140 
AL(14 
400)  , 
400)  , 
GER  P 
NSION 

1=1 

K=1 
0   J  =  1 

(5,1 
E  (6, 
INUE 

READ 

WRIT 

IF     ( 

INUE 
1=1  + 

INUE 
READ 
WRIT 


IN    WEARMETAL 


ANALYSIS******* 
2ES    THE    10    HOUR    TRENDS 
UER 

R(1400)  ,  DEL  FE(  1400)  ,  DELAL  (140  0) 
0)  rDELCU  (1400)  .DELPB  (1400)  . 

0)  , TENDCR(  1400)  ,  TRNDPB  ]l4  0l 


E(HQQ)  *AL  .(1400).,  CR  .(1400)  , 


0) 


1400)  ,TRNDCU(1400)  ,  TRNDFE  ( 1400) 


CU 
RE 
NAME(75) 


,75 

10)  NAME  (J) 
111)  NAME(J) 


(5,  100)  HR  (I)  ,FE  (I)  ,  AL  (I)  ,CR  (I)  ,CU  (I)  ,PB  (I) 
E  (6,200)  HR(I)  ,FE(I)  ,AL<I)  ,CR(t)  ,CU  (I)  , 

HR(I)  .EQ.  (0  .)  )  GO  TO  5 


1 


(5,  100)  HR  (I)  ,FE  (I)  ,  AL  (I)  ,CR  (I)  ,CU  (I)  ,FB  (I) 
E  (6,200)  HR(I)  ,FE(I)  ,AL  (I)  ,CR(I)  ,CU(I)  , 


WRITE  (6, 
*EELCU 

CONT 
DELH 
DELF 
DEL  A 
DELC 
DELC 
DELP 
WRIT 
*DELCU  (K) 

continue 

tend: 

TEND 

TFNDi 

TEND' 

trnd: 

writ: 

*tfnd 

K  =  K  + 

GC  T' 

ccnt: 

1=1- 

cont: 

CALL 

writ: 

CALL 

writ: 

CALL 

writ: 

CALL 
WRIT 


DELAL  (K)  ,DELCR(K)  , 


DELP 


DELHR  (K)  ,DELF 


,  DELAL  (K)  ,  DELCR(K) 


DEIFE 
DELAL 
DELCR 
DELCU 
DELPB 
TEND] 


K)  /DELHR 
K)  /DELHR 
K'i  /DELHR 
K  /DELHR 
K)  /DELHR 


K) 

3 

K) 


*10, 
*10 

*10 

*10. 
*1Q 


E(K)  ,TRNDAL(K)  ,TRNDCR  (K)  ,TRNDCU  (K) 


GP     (TRNDFE, K,0) 

210) 

GP     (TRNDAL,K,0) 

220) 

GP     (TBNDCR,K,0) 

230) 

GP     (TRNDCU,K,0) 

24  0) 
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CAIL  KESTGP  (TENDPB,  K,0) 

WRITE  (6,250) 
C  HBITE  (6,210) 
C      WHITE  (11.*)  TENDFS 


C  WRITE  (6,220) 

C  WBITE  (7,*)  TENDAL 

C  WRITE  (6,230) 

C  WHITE  (8,*)  TENDCR 

C  WHITE  (6,240) 

C  WRITE  (9,*)  TRNDCU 

C  WHITE  (6.250) 

C  WHITE  (10,*)  IRNDPB 
STCP 

100  FORMAT  (1X,6F10.5) 

110  FORMAT  (1X,A10) 

111  FORMAT  (1X,A10[ 
C200  FCEMAT  (1X,6F10.2) 

210  FORMAT  ( 1 X, • HISTOGRA MS  OF  PIECEWISE  TRENDS  FOR 


*IECNf  ) 
FORMAT 
*ALCMINUI 


220   FORMAT  ( 1 X, • HISTOGRA MS  OF  PIECEWISE  TRENDS  FOR 

~IH"  ) 
230   FORMAT  ( 1 X, • HISTOGRA MS  OF  PIECEWISE  TRENDS  FOR 

♦COPPER1) 
240   FORMAT  ( 1 X, • HISTOGRAMS  OF  PIECEWISE  TRENDS  FOR 

♦CHROMIUM1 ) 
250   FORMAT  ( 1 X, ■ HISTOGRA MS  OF  PIECEWISE  TRENDS  FOR 
*LEAD«) 
END 
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c  *****  ERCGRAM  ANALGEN  ** 
C  PROGRAM  TO  GENERATE  SIMU 
C  RICHABE  F.  BAUER 

REAL    HR  (600)  ,  IRON  (60 

*DELIRN(599)  ,TRDIRN(5 

*DFLIE1j2396)  ,TEDIRT( 

INTEGER     PEErPEEI,PRE 

SEFD=10. 

CALL    LNORM(SEED,ERRO 
DC    5    1=1,550 

HR  (I)=10.*(FLOA 
IRON  (I)  =  (  (.3648 
J=I  +  600 
K=I+1200 
L=I-H800 
IRON1  (I 
IRON1  (J 
IRON1  (K 
IRON1  (L 
5       CONTINUE 
B0=IRON (I) 
SOI=IRON1  " 
BCJ=IRON1 
E0K=IRON1 
BOL=IRON1 
DC    7    1=55 


10 


C200 


Jj*  *fr  Jfr    j«^   -*jfc 

LATED  DATA  FOR  USE  IN  ANALYSIS 

0)  -ERROR  I  2  40  0)  ,  DEL  HR  (5  9  9)  , 

99)  .IRON!  (2400)  , 

2396) 

J,PREK,PREL 

R, 2400, 2,0) 

T(I)) 

) *HR  (I))  +ERROR (I) 


*HR  (I 
*HR  (I 
*HR  (i; 

*hr  \r 


♦ERROR  (I) 
♦ERROR  (J 
+  ERROR  (K 
♦ERROR  (L 


600 


jj=i-5oo 

HR  (I)=10.*(FLOA 

J=I  +  60  0 

K=I+1200 

L=I+1800 

IRON  (I)  =  E0+(  (.7 

IRONT(I)  =  E0I  + 

IRON1  (J)  =  E0J  + 

IRON1(K'i  =  E0K  + 

IRON1  (L    =E0L  + 
CONTINUE 
DO    1C   JJ=2,600 

I=JJ 

J=JJ+600 

K=JJ+1200 

L=JJ  +  1800 

PRE  =  JJ-1 

DELHR(JJ)=HR  (JJ 

DELIRN  (JJ)=IRON 

TRDIRN  (JJ)=(DEL 

PREI=I-1 

PRSJ=J-1 

PREK=K-1 

PREL=L-1 

DELIR1 

DELIR1 

DELIR1 

DELIR1 

TRDIE1 

TRDIR1 

TRDIR1 

TRDIE1 

WRITE 
CONTINUE 
CALL    HISTGP 
WRITE    (4,*) 
CALL    HISTGP 
STCP 

FORMAT  (1X,1F10.5) 
END 


T(JJ)) 


296)*HR(I)  ) +ERROR  (I) 
.7296)  *HR  (I)  )  +  ERROR(I) 
.9554)*HR  (I)  TERROR  (J 
1.0282)  *HR  (if  )  +ERROw  (K) 
.66  12) *HR  (I) )  +ERROR(L) 


)  -HR  (PRE) 

(JJ)  -IRON  (PRE) 

IRN  (JJ)/DELHR  (JJ)  )  *10 


(I 
J 

K 
L 

ij: 

K)  = 

,L    = 

6,2 


=IRON1 

=IRON1 

=IRON1 

=IRON1 

DELI 

DELI 

DELI 

DELI 


50)     T 


IRON1  (PREI) 
IRON1 (PREJ 
IRON1 (PEEK) 
IRON1 (PREL) 
/DELHR  (JJ 
/DELHR  (JJ 
/DELHR  (JJ 
/DELHR  (JJJ 
N  ( JJ) 


*10. 
*10. 

♦  10. 

♦  10. 


(TRDIRN 
TEDIRN 
(TEDIR1 


PRE, 15) 
2396,15) 
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